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Abstract 

Lattice model with long-range interaction of power-law type that is connected with 
difference of non-integer order is suggested. The continuous limit maps the equations of 
motion of lattice particles into continuum equations with fractional Griinwald-Letnikov- 
Riesz derivatives. The suggested continuum equations describe fractional generalizations 
of the gradient and integral elasticity. The proposed type of long-range interaction allows 
us to have united approach to describe of lattice models for the fractional gradient and 
fractional integral elasticity. Additional important advantages of this approach are the 
following: (1) It is possible to use this model of long-range interaction in numerical simu¬ 
lations since this type of interactions and the Griinwald-Letnikov derivatives are defined 
by generalized finite difference; (2) The suggested model of long-range interaction leads to 
an equation containing the sum of the Griinwald-Letnikov derivatives, which is equal the 
Riesz’s derivative. This fact allows us to get particular analytical solutions of fractional 
elasticity equations. 

PACS: 45.10.Hj; 61.50.Ah; 62.20.Dc 


1 Introduction 

Discrete system of long-range interacting particles serve as a model for numerous applications 
in mechanics and physics m ia si a ia ei in si ei mu hb. Long-range interactions are important 
type of interactions for complex media with non-local properties (see references in [12]). 

Using the fractional calculus pun], we consider long-range interaction of a special type 
to describe a fractional generalization of elasticity theory. We transform the set of equations 
of motion of coupled individual lattice particles into the equation of non-local continuum that 
contains spatial derivatives of non-integer orders. It allows us to consider different lattice models 
for generalization of elasticity theory by applying methods of fractional calculus. The theory of 
fractional derivatives and integrals has wide applications [15, 16j 17. Ml [El 20. EU122, 12] and 
it is a powerful tool for the analysis of different nonlocal continuum models with nonlocality of 
power-law type. 
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Non-local continuum mechanics has been treated with two different approaches: the gradient 
models (weak non-locality) and the integral models (strong non-locality). The correspondent 
constitutive relations have the form 

&ij C'ijkl^kl A l^kl) , 

where C t] ki is the elasticity tensor, and a %] are the strain and stress tensors, respectively. The 
additional parameter l s is an internal length scale that can be defined by lattice constant. The 
operator C is the integral or differential operator of integer order [31[ |32, 3U [33], 35] • Recently, 
the problems of non-local continuum are described by means of fractional calculus. Fractional 
models of integral non-local elasticity are considered in [23 25] 23 121,23 221 3D] - In fractional 
elasticity models, the operator C is the integral or differential operator of non-integer order. To 
generalize continuum equations by using fractional calculus, we should represent these equation 
through the dimensionless coordinate variables. Therefore the coordinates x, y, z, the vector 
r, the value r = |r|, and the parameter l 2 are dimensionless in the suggested fractional models. 

In this paper we consider one-dimensional lattice model with long-range interactions of 
of Griinwald-Letnikov-Riesz type [12] • A feature of suggested long-range interactions is that 
the interactions terms have the form of fractional differences of non-integer orders. The first 
advantage of this approach is based on the properties of the Griinwald-Letnikov fractional 
derivative [13 23. The suggested type of long-range interaction of lattice particles allows us 
to have united description of lattice models for the fractional gradient and fractional integral 
elasticity that is characterized by the non-locality of power-law type. The second important 
advantage of suggested approach is the ability to directly use the model of long-range interaction 
in numerical simulations since this type of interactions and the Griinwald-Letnikov derivatives 
are defined by generalized difference of non-integer order. We assume that the suggested form 
of long-range interaction can be used for different scheme of simulations in fractional gradient 
and integral elasticity models. 

There are some problems with application of the Griinwald-Letnikov derivatives, since there 
are currently very limited number of analytical solutions for differential equations with frac¬ 
tional Griinwald-Letnikov derivatives in contrast to the equations with derivatives of Riemann- 
Liouvillc, Caputo and Riesz types [IT]. We note that the suggested model of long-range inter¬ 
action leads to an equation containing the sum of the Griinwald-Letnikov derivatives, which is 
equal the Riesz derivative. Using this connection, we derive some particular analytical solutions 
p3] for fractional elasticity equations with the Riesz derivatives of non-integer orders. 

2 Griinwald-Letnikov fractional derivatives and integrals 

Fractional differences of non-integer orders 

The Griinwald-Letnikov derivatives have been introduced by Griinwald in 1867 and indepen¬ 
dently by Letnikov in 1868. Definition of the Griinwald-Letnikov fractional derivatives are 
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based on a generalization of the usual differentiation of a function f(x) of integer order n of 
the forms 


D n J(x) 

D n J(x) 


lim 

o 


lim 

o 


A hf(x) 

h n 

h n 


J 


J 


where AJJ and are forward and backward finite differences of order n of a function f(x) 
with a step h and centered at the point x. The nth-order forward and backward differences are 
respectively given by 


xnx) = Y,i-i) t 

k =0 


f{x+(n-k)h), 


( 1 ) 


v;;/m = X>i) 1 (?)/(*-**)• (2) 

k =0 ' ' 

Note that 

A n h f{x) = (-lJTJW. (3) 

The difference of a fractional order a > 0 is defined by the infinite series (see Section 20 in [13j ) 


OO / \ 

= ( 4 ) 

where the binomial coefficients are 

/ a\ T(a + 1) 

U/ = T{k + l)T(a-k + iy 

For h > 0, the difference (j4j) is called left-sided fractional difference, and for h < 0 it is called 
a right-sided fractional difference. We note that the series in ([4j) converges absolutely and 
uniformly for every bounded function f(x ) and a > 0. 

For the fractional difference, the semigroup property 

V2Vf/M = (OO, ()>0) (5) 


is valid for any bounded function f(x) (see Property 2.29 in [14] page 121). 
The Fourier transform of the fractional difference is given by 


T{X7 a h f(x)}(k) = (1 - exp {ikh}) a F{f(x)}(k) 


for any function f(x) G Li(M) (see Property in [H] page 121). 
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Grunwald-Letnikov fractional derivatives 

The definitions (P) and (P are used to define the Grunwald-Letnikov fractional derivatives by 
replacing n G N in by a > 0. The value h n is replaced by h a , while the finite difference is 
replaced by the difference of a fractional order a. 

The left- and right-sided Griinwald-Letnikov derivatives of order a > 0 are defined by 


GL D a x± f(x) 


lim 

h —^0“b 


V±hf( x ) 

\h\ a 


Note that the Grunwald-Letnikov derivatives for integer orders a = n G N is 


( 6 ) 


gl di ± Hx) = (-iyd\ 


p) 


We also note that these derivatives coincide with the Marchaud fractional derivatives of order 
a > 0 for f[x) G L p (M), 1 ^ p < oo (see Theorem 20.4 in [13]). The properties of the 
Griinwald-Letnikov fractional derivatives is described in Section 20 of the book [13]. Then (}6j) 
can be represented by the Marchaud fractional derivatives 


GL 


D a x± f(x) 


a 

T(1 - a) 


[°° f(x)-f(x T z) 

Jo - Q+1 


if f(x) G L p (M), where 1 < p < 1/a and 0 < a < 1. 


Griinwald-Letnikov fractional integral 

It is interesting that series (J4|) can be used for a < 0 (see Section 20 in H) and equation (J6]) 
defines the Grunwald-Letnikov fractional integral if 

\f(x)\ < c(l + \x\)~^, fi> |a|. (8) 

The existence of the Grunwald-Letnikov fractional integral means that we have a united def¬ 
inition of fractional derivatives and integrals. It allows us to have a united approach to the 
fractional gradient and integral elasticity. 


Commutativity and associativity of Grunwald-Letnikov derivatives. 

In [23J has been showed that, for analytic functions, the fractional derivatives of the Griinwald- 
Letnikov type have some nice and useful properties and a semi-group structure. It is connected 
with the fact that differences of fractional order satisfy the semigroup property (J5]). Using the 
semi-group property 

GL D a x± gl D/ ± = GL D«t\ (a > 0, (3 > 0). (9) 
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Ortigueira, Rivero, Trujillo proved [23] that the Griinwald-Letnikov derivatives are commutative 
and associative operators. These properties are very important to application in mechanics, 
but only a few types of fractional derivatives have them. These properties allow us to represent 
a single fractional derivative, which appears in the equation of the lattice model, as a product 
of derivatives, and thus get the fractional gradient and integral elastic constitutive relations. 


Grunwald-Letnikov-Riesz fractional derivative 


We can define a fractional derivative of order a > 0 by 


aLR D“f(x) 


i v;/(x) + VVW 

-7-TT hill -;-;- 

2cos(a;7r/2) \h\ a 


( 10 ) 


This derivative coincide (see Section 20.1 in [13] ) with the Riesz fractional derivative of 
order a > 0 if a ^ 1, 2, 3,..., i.e. we have the relation 

GLR D a j(x) = R n a j(x ), (11) 


where 


R 


D a J(x) = 


a 


f(x + z ) -2 f(x) + f(x-z) 


2T(1 — a ) cos(ck7t/2) 


yOi~\~ 1 


dz. 


( 12 ) 


Therefore the fractional derivative is called the Grunwald-Letnikov-Riesz derivatives of 
order a >0 (a^ 1,2,3,...) in [ IB] . 


3 Lattice with linear nearest-neighbor interaction 


Let us consider equations of motion for particles with the nearest-neighbor interaction of the 
form 

r \2 

M —u n (t ) = KA k n u k (t) + F n (t), (13) 

where 

771=1 

Ku m {t) = ^(-l) m+1 (u n+m {t) +« n _ m (i)). 

m= 0 


For these equations we can give the well-known statement regarding the nearest-neighbor in¬ 
teraction (see for example [3&, 371 18]). 

Proposition 1. In the continuous limit the lattice equations of motion 11 dl) maps into the 
continuum equation 


d 2 u(x,t ) _ 2 d 2 u(x,t ) 1 

dt 2 e dx 2 p 


(14) 
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with the mass density p = M/Ah, the Young’s modulus E = Kh/A, the force density f(x,t ) = 
F(x,t)/Ah, the cross-section area of the medium A and the inter-particle distance h, where 


r,2 E Kh?_ 

e p M 


(15) 


is a finite parameter. 

Proof. To derive the equation for the field u(k, t ), we multiply equation (IT3|) by exp(— iknAx), 
where Ax = h and summing over n from — oo to +oo. Then 

+OO rj2 + OO + OO 

M Y e ~ iknAx -^r = K e- iknAx [u n+1 -2u n + u n . 1 ]+ Y e~ iknAx F(n,t). (16) 

n =—oo n =—oo n =—oo 


The first term on the right-hand side of (TT6T) is 


e ~ lknAX (wn+l - 2u n + U n _i) = 


= e 


ikAx 


-(-oo 

E 

l= — C 
+C 

E 


+oo 


+oo 

e~ iknAx u n+1 - 2 Y 

n =—oo 
+oo 

e -ikmAx 


e~ iknAx u r 


+oo 

E 


e~ iknAx u n _x = 


+oo 


2 Y e~ iknAx u n + e~ ikAx Y 

m =—oo n =—oo s =—oo 

= e ikAx u(k, t ) - 2 u(k, t ) + e~ ikAx u(k , t) = 

= [e ikAx + e~ ikAx - 2 ]u{k,t) = 2^cos(A:Ax) - l)u(k,t) = -4 sin 2 
Here we use the definition of u(k,t ) on [—/2, fc 0 /2] by the equation 

+oo 

t(M) = Y u nif) e~ lkXn = T A {u n (t)}, 


e~iks Ax u _ 


(¥) 


u(k, t). 


U( 


(17) 


and x n = n h, where h = 27t//c 0 is distance between equilibrium positions of the lattice particles. 
As a result, we have 


M d ~ ui : k : ^ = K J a (kAx) u{k , t) + F A {F n {t)}, 


dt 2 


where 


J a (kAx) = —4 sin 


2 i kAx\ 


2 J 


(18) 

(19) 
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For h = Ax —y 0, the asymptotic behavior of the sine is sin(fcAx/2) ~ kAx/2. Then (1T9|) can 
be represented by 

J a {kAx) ph — ( kAx ) 2 . 

Using the hnite parameter C 2 = K h 2 /M, after division by the cross-section area of the medium 
A and the inter-particle distance h, the transition to the limit h = Ax —> 0 in equation (fT8l) 
gives 

AA *> = _ C * eu(k, t) + t)}, ( 20 ) 

where we use 0 < | C 2 \ < oo. The inverse Fourier transform J 7-1 of (1201) has the form 


d 2 J 7 ~ 1 {u(k,t)} 
dt 2 


- C 2 T 1 {k 2 u(k,t)} + ^ 


Then we can use the connection between second derivative and its Fourier transform in the 
form k 2 <—» — d 2 jdx 2 . As a result, we obtain the continuum equation (TT4|) . □ 


As a result, we prove that equations (fT5]i give the continuum equation with derivatives of 
second order only. 


4 Lattice with Grunwald-Letnikov-Riesz long-range in¬ 
teraction 

In this section we describe the type of long-range interaction that is suggested in [12] (see 
Section 8.19). Let us consider a lattice system of interacting particles, whose displacements 
from the equilibrium are u n (t), where n £ Z. We assume that the system is described by the 
equations of motion 

d 2 

M d^ Un ^ + 9 B n( a ) U m(t) - F n (t ) = 0, (21) 

where g is the coupling constant for long-range interaction that have the form 

+oo 

B™{a)u m {t) = ^2 ba(.m) ( u n+m {t) + u n - m (t)j , ( 22 ) 

m =0 


and the function b a (m) is 


(-l) m 

b a (m) = — - ——-— T- 

1 (m + 1) f (a — m + 1) 


(23) 


This type of long-range interaction has been called the Grunwald-Letnikov-Riesz interaction 
(see Section 8.19 in [12]). Let us give the main statement regarding this interaction. 
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Proposition 2. In the limit h —> 0 lattice equations pi]) . ( OH]) rozf/i / f22]) give the continuum 
equation 


d u{x, t) + c{a) GLR Df u(x, t ) - -f(x, t) — 0, (a6R,a/ ±1, ±3, ±5, 

at" p 


(24) 


where 


is a finite parameter, and 


C(a ) = 


2 cos(ct 7r/2) g /i° 
T(a + 1)M 


(25) 


GLR r^a 


D a = - 

x 2 cos(a7r/2) 




(26) 


is t/ie Griinwald-Letnikov-Riesz fractional derivative of order a, andu(x,t ) is a smooth function 
such that u(nh,t ) = u n (t). 

Proof. We define smooth functions u(x,t ) and F(x,t) such that 


u(nh,t ) = u n (t), F{nh,t ) = F n {t). 
Then equations (]2T]h (1221) with (1231) can be represented as 


d 2 u(x,t) gh‘ 


dt 2 


+ — > b a (m) -— (u(x + mh, t ) + u(x — mh , t) ) — —Fix, t ) = 0. 
M ^ \ 1 ' v V M v 


(27) 


m=0 


After division by the cross-section area of the medium A and the inter-particle distance h it is 
found that 


d 2 u(x,t) C(a ) 


<9t 2 


+oo 

£ 


(—l) m r(cK T 1) w(a: + mh, t) + u(x — mh, t) 1 


2 cos(ck7t/2) ^ r(a; — m + 1) T(m + 1) 


h° 


P 


f(x,t) = 0 
(28) 

with the mass density p = M/Ah, the force density f(x,t ) = F(x,t)/Ah, and C(a) is defined 
by (1251) . Using the definitions of the left-sided and right-sided fractional differences and the 
limit h —> 0+, we obtain 


d 2 u(x, t ) C(a) 


Vfu(x,t) + Vf h u(x,t) 1 n 

„ „ , . , , hm -—- tlx, t) — (J. 

dt 2 2 cos(a7r/2) o \h\ a p 


(29) 


Using the Griinwald-Letnikov-Riesz derivative (ITU)) , equation (122]) can be rewritten in the form 












5 Fractional gradient and integral elasticity of Grunwald- 
Letnikov-Riesz type 


Let us consider a system of interacting particles, whose displacements from the equilibrium are 
u n {t), where n £ Z. We assume that the system is described by the equations of motion 


0 2 U n {t) 

d t 2 


- KA£u m {t) + gBZ{a)u m {t) - F n (t) = 0. 


M 


In the limit h —> 0 equations (I30T) gives the continuous medium equation 


d 2 uix,t) 2 d 2 uix,t) 


- Cl + C(a) GLR D X u{x , t) - - fix, t) = 0, 


dt 2 e dx 2 

where Cl is defined by (IT5lh C(a) is defined by (l25]h and 


1 

P 


(30) 


(31) 


a £ K, =tl, ±3, ±5,... . 

Let us consider equation f|3Tjl for two cases: a > 2 and 0 < a < 2. 

For a > 2, we can use gl D]. ± = ±Ld, and the semi-group property (J9]) , to represent the 
Griinwald-Letnikov derivatives in the form 


GL 


D a 

U x± 


= D 


1 GL D«± 2 Dl 


ia > 2 ). 


(32) 


Therefore, we have 


GL 


D a x+ + gl D x _ = Dl 


GL 

u x+ 


GL jdd —2 


+ 


Dl 


(33) 


As a result, we can use 


GLRr^a _ 7~\1 GLR T^oi —2 1 

-* — ' rv» J- — ' rv* J — ' rp -*■ —^ 


in order to rewrite equation d3T]) in the form 


(34) 


—~h ^ - c l £> 2 «(^.«) + C(a) G “r>r 2 = (35) 

The correspondent constitutive relation for the continuum equation ()35l) can be derived by the 
momentum balance equation 


P 


9 U ^2 — = D x a ( X ) + /OM), 


and the strain-displacement relation for small deformations 

e(x) = D\u(x). 


(36) 


(37) 
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The commutative and associative properties of fractional Griinwald-Letnikov derivative and 
relation (j32l) allow us to represent a single fractional derivative in equation (l24lh as a product 
of derivatives, and thus get the fractional gradient constitutive relations. 

Using (136|) and (J37]), equation (135|) gives the fractional constitutive relation in the following 
form 

cr(x,t ) = E (e(x,t) =F l 2 (a) R Y) X ~ 2 e(x,t)) ( a > 2), (38) 

where E is the Young’s modulus, is the Riesz fractional derivative of order a — 2, which 

is equivalent by (HIT) to the Griinwald-Letnikov-Riesz fractional derivative GLFt D " -2 . As a 
result, we can state that equation (I3D with a > 2 describes a fractional gradient elasticity. In 
equation (|38|) . we use 

\C{a)\p 2 \g\ li a ~ 2 | cos(a7r/2)| 


ls ^ E K T(a + 1) 

is the scale parameter of fractional elasticity. Sign in front of this scale parameter Z 2 ( 


(39) 
a) in 

equations (138|) is determined by the sign of the expression g cos(a7r/2), i.e. the sign of the 
coupling constant g of lattice vibrations and the value of the order a of long-rang interactions. 
If g cos(a7r/2) > 0, then we get the minus in front of l 2 (a ) in (l39lh For the case a = 4 and 
g > 0, we derive the constitutive relations for the gradient elasticity model with the minus 
in front of l 2 . We get the equations for phenomenological gradient elasticity model from the 
lattice model equations. Note that normally it is considered a phenomenological model does 
not have a corresponding microscopic model [35]. 

For the case 0 < a < 2, we cannot use the properties (1321) and (|34li because semi-group 
relation (J9]) holds for positive orders. For 0 < a < 2, we can use the relation (11111 and the 
Riesz’s analytic continuation (see equation (17) in [46]) of fractional integrals of order a to 
negative values of orders a > —p in the form 


iija _ ( — R\a+2p 


(40) 


where pGN and R \ a+2p is the Riesz fractional integral of order a + 2 p. Using equations (17), 
(21) and (22) from the Riesz’s review paper [46], we get 


*D“ 


-Rj2 —a 


A = A R I 4 ~ a A. 


For one-dimensional case, we have A = (D 4 ) 2 , and we can use the relation 

R-r^oi _ A -Rt4— a A _ nl t2-o nl 


= a K r~ a a = d l x r~ a d‘, 


(41) 


(42) 


where J 2 a is the fractional integral operator of order (2 — a) that is defined by 


r2—a 


= Di T 


1 jj 1 


(43) 


As a result, equation (l3l| with a < 2 describes a fractional integral elasticity. Using 
(j37D and (142|) . equation (13TT) gives the fractional constitutive relation in the following form 


<7 


(x, t) = E ( e(x , t ) =F l 2 s {a) R J 2 X a e(x, t )) (0 < a < 2), 


(44) 
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where J 2 ~ a is the fractional integral (l43li of order 0 < 2 — a < 2, and l 2 s {a) is defined by (1391) . 
Sign in front of this scale parameter l 2 (a) in equations (j44h is also determined by the sign of 
the expression g cos(a7r/2). 

Let us give a remark about the scale parameter l s (a). Equation (139|) can lead to incorrect 
conclusion about the behavior of the parameter l 2 (a ) for h —> 0 in the case 0 < a < 2. Using 
(USD, the dimensionless parameter (15§P can be written as 


2 \g\ h a | cos(a7r/2)| 


Cl M T(a + 1) 


(45) 


Because the value of C 2 if hnite, then behavior of the parameter l 2 (a) for h —> 0 has the 
identical type for a > 2 and 0 < a < 2, such that l 2 (a) is proportion to h a . 

As a result, equations (|38|1 and (1441) describe the fractional gradient and fractional integral 
elasticity of non-local continuum. If 0 < a < 2, we have a fractional integral elasticity, and if 
a > 2, then equation describes fractional gradient elasticity. It can call the fractional elasticity 
of Griinwald-Letnikov-Riesz type. 


6 Solutions of fractional elasticity equations 


Let us consider more general lattice system of interacting particles that is described by the 
equations of motion 

M + K Ku m {t) + Y,g k Bn^k)u m {t) = F n (t). (46) 

fc=i 

The correspondent continuum equation for the fractional elasticity of the Griinwald-Letnikov- 
Riesz type has the form 

- U fa2 ^ ~ C e D l U ( X > t )+'%2 C ( a k) R V* k u{x) = -f{x). (47) 

fc=l p 


Here to get analytical solution of the fractional equations of nonlocal elasticity, we use the 
equivalence (see Section 20.1 in [13]) of the Griinwald-Letnikov-Riesz derivative GLB D% and 
the Riesz derivative in the form 


GLR 


DZ = 


2 cos(cK7r/2) 


GL 


D a x+ + 


GL 


DZ ) = 


(48) 


Equation (l47ll can be considered as a fractional generalization of the higher-order strain-gradient 
models [34J and the higher-order integral elasticity models. In the static case ( D 2 u(x,t ) = 0) 
the fractional elasticity equations (|47li is 


E C (“*) = -/W, (49) 

P 

t— 1 ' 
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where x £ R, N G N and m ^ 1, C(ptk) G R, cun > ....an > 0 

Equations (j47T) and (j49li involve the one-dimensional Riesz fractional derivatives given by 


r d xn 


i r (A»tn 

di(l,a t ) |z|“» +1 iiz ’ 


(50) 


where a& < /, k = 1,2, ...,m, and (A ™f)(z) is a finite difference of order m of a function /(a:) 
with a vector step zeR and centered at the point i£l: 


(V"/)M = El- 1 )' 


m\ 


k =o 


k\(m — fc)! 




fcz). 


The constant d n (m, a) is defined by 


, , . 7r 3 / 2 A m (a) 

l[m ’ “ 2°T(1 + a/2)T(l/2 + a/2) sin(7ra/2) ’ 


where 


P) j!(m-j)! 


Note that the hypersingular integral D“/(x) does not depend on the choice of m > a. 

Equations (|49|) are solvable, and it particular solutions are given (see Theorem 5.25 in [1%] ) 
by the formula 


1 r+oo 

u(x) = - G a (x - z) f(z)dz, (51) 

P J —oo 


where 


G r 


x 


Z = - 

IT 


J2C(a k )\ ak + C*\ 2 

^ k =1 


cos(A|x|) d\. 


(52) 


Let us consider one-dimensional W. Thomson (1848) problem [39]. We determine the de¬ 
formation of an infinite elastic continuum, when a force is applied to a small region in it. We 
consider one-dimensional elastic media with power-law nonlocality that is described by the 
equation 

- Cl Dlu(x) + C(a) h D;«(x) = - f(x). (53) 

P 

Note that equation (}53]) coincides with equation (13T]) for static case. If we consider the defor¬ 
mation at positions x, which are larger compare with the size of the region, where the force is 
applied, we can suppose that the force is applied at a point. In this case, we have 


f(x) = f 0 5(x). 


(54) 
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Then the deformation, which is a particular solution of equation (|53lh will be described by the 
equation 


u{x) 


fo_ f°° cos(A|x|) 
irpJo C% A 2 + C(ot)\ a ' 


(55) 


If 0 < a < 2, the solution (T55|) corresponds to the fractional integral elasticity, and if a > 2, 
then the solution (T55|) corresponds to the fractional gradient elasticity. 

We can consider more general model of lattice with long-range interaction in R 3 , where all 
particles are displaced from its equilibrium in one direction, and the displacement of particles 
is described by a scalar field u(x), where x = |r| and r e R 3 . The correspondent continuum 
equation of the fractional elasticity model is 


- Au(x) + ^((—A Yl\)(x) = T-./M, 


(56) 


where (—A) a//2 is the fractional Laplacian pH], The displacement vector u(x) of the point force 
(|54|) has the following form 


u(x) 


1 fo 

47rU 2 p x 


C 2 


(57) 


where x — |r|, and 


C 2 , a (x) 


2 r°° A sin(Ax) 
nj 0 A 2 + ( C{a)/Cl)\ 


(58) 


We note that the asymptotic behavior x = |r| —» 0 of the scalar field u{x) does not depend 
on the parameter a. Using (see equation (1) of Section 2.3 in the book [38]), we obtain the 
asymptotic behavior (x —> oo) for C 2: a (x) with a < 2 in the form 


a 


2 ,a 


» = - 

VT 


A sin (Ax) 


A 2 + (C(a)/C*)\< 


dX ~ A Q (a)- 


x 


2—a 


oo 

E 

k =1 


A 


k(a 


x (2-a)(k+l) ’ 


(59) 


where 


A o(®) = 2 ^ e T(2 - a) sin , (60) 

7TC{a) \2 / 

o / f^2 \ k+1 poo 

M a ) = — ( 7 7 TT) / zC- 01 ^- 1 sm(z)dz. (61) 

71- \C (OL) J J o 

As a result, we have in the framework the fractional elasticity, the displacement field of the 
point force in the infinite media with non-locality of power-law type is given by 

“<*> " (0<Q<2) (62) 


for the long distance x ^ 1. 
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In Figures 1-3 we present some plots of the factor C 2 ,o(x) = exp(— x) and factors C' 2 , a {x) 
with C{a)/C g = 1 for different orders of 0 < cx < 2, i.e. for the fractional integral elasticity. 
The values of the factors C 2 , a (x) and C 2 ,o(x) are plotted along the Y-axis, and the values of 
the position x — |r| are plotted along the X-axis. 

For the case cx > 2, i.e. for the fractional gradient elasticity, the asymptotic behavior 
x = |r| —» oo of u(x) does not depend on the parameter cx. The asymptotic behavior of the 
displacement field u(|r|) for |r| —>■ 0 is given by 


u{x) ~ 


/o T((3 — cx)/2) 

2" 7 T 2 yfH p C (a) r(a/2) 


(2 < cx < 3), 


(63) 


U(X 


fo 


(cx > 3). 


(64) 


2-kq>p(CI) 1 3 /“ C 3 /“(cx) sin(37r/cx) 

Note that the function C- 2 t a for the fractional gradient case (cx > 2) has a maximum. If ex = 4 
then we have the well-known case of gradient elasticity [ 35] . 

In Figures 4-6 we present some plots of the factor C 2 ,o(x) = exp(— x) and factors C 2 , a (x) 
with C[ol)/C 1 = 1 for different orders of 2 < cx < 6, i.e. for the fractional gradient elasticity. 
The values of the factors C 2 lCe {x) and C 2 ,o(x) are plotted along the Y-axis, and the values of 
the position x — |r| are plotted along the X-axis. 


7 Conclusion 

A lattice model with long-range interaction of Griinwald-Letnikov-Riesz type is suggested. In 
the continuum limit we derive continuum equations with spatial derivatives of non-integer 
order cx. The correspondent constitutive relations describe fractional generalization of gradient 
elasticity for cx > 2 and fractional integral elasticity for 0 < cx < 2. The suggested lattice model 
is considered as a microscopic model of the fractional non-local elastic continuum. We can 
note that a fractional nonlocal continuum model can be obtained from different microscopic 
or lattice models. The benefits of suggested formulation of fractional elasticity are following. 
Firstly, the Griinwald-Letnikov-Riesz derivatives in the fractional continuum equations are 
defined by fractional differences. It can be directly used in numerical simulations of fractional 
gradient and fractional integral elasticity models. Secondly, the suggested type of long-range 
interaction for lattice particles allows us to have united lattice model for the fractional gradient 
and fractional integral elasticity. We assume that the suggested approach can be generalized 
for 3-dimensional case, for finite strains and plasticity. An extension of the suggested model for 
these cases can be realized by the methods suggested in [12] (see Sections 8.2, 8.14, 8.15) with 
some modifications. In addition, we note that the model of fractional gradient and integral 
clastic continuum has an analog in the plasma-like dielectric material with power-law spatial 
dispersion [4DJ SI]. Fractional models of complex material with power-law non-locality allows 
us to predict unusual properties of materials that are characterized by long-range inter-particle 
interactions. These materials can demonstrate a common or universal behavior in space by 
analogy with the universal behavior of low-loss dielectrics in time [42j SB] SH S5]. 
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th a = 0.7 is drawn with black color and the 
color, where x — |r| and C(a)/C | = 1. 



th a — 1.1 is drawn with black color and the 
color, where x — |r| and C(a)/Cl = 1. 






,(x) with a = 2.2 is drawn with black color and the 
y gray color, where x — |r| and C(a)/Cl = 1. 
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Figure 5: Plot of the function y = C 2 , a { x ) with a = 3.6 is drawn with black color and the 
factor y = C 2 ,o(^) = exp(— x) is drawn by gray color, where x — |r| and C{oi)/C\ = 1. 



Figure 6: Plot of the function y = C < 2 ,a(^) with a = 5.9 is drawn with black color and the 
factor y = C 2 fi(x) = exp(— x) is drawn by gray color, where x — |r| and C{a)/C\ = 1. 
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